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ABSTRACT 

An  examination  is  made  of  the  physical  bases  for  the  description 
of  turbulence  by  probability  distributions,  and  for  the  closure  problem 
which  arises  from  this  description.   Then  a  discussion  is  given  of  the 
analytical  approaches  to  the  closure  problem  which  have  been  proposed 
by  Deissler,  Hopf,  Kraichnan,  Proudman  and  Reid,  and  Tatsumi .   Some 
mathematical  inter-relations  of  these  schemes  are  cited.   Certain  con- 
sequences of  the  several  approaches  are  illustrated  by  applying  them  to 
an  idealized  turbulent  convection  problem  for  which  they  all  yield  ex- 
plicit solutions. 
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1.   INTRODUCTION 

Turbulent  flows  vhlch  occur  in  nature,  or  the  laboratory,  can  appear 
Infinitely  complex  to  the  eye.   This  is  faithfully  reflected  in  turbulence 
theory.   A  central  feature  of  the  theory  is  an  infinite  set  of  dynamical 
equations  which  couple  together  all  the  moments  of  a  statistically  distri- 
buted velocity  field.   The  closure  problem  associated  with  these  equations 
may  be  posed  as  follows:   How  are  they  to  be  replaced  by  a  finite  set  that 
yields  the  limited  statistical  information  of  actual  interest?  This  ques- 
tion must  be  answered,  explicitly  or  implicitly,  in  order  to  solve  most 
problems  in  turbulence  theory.   Consequently,  the  closure  problem  has  re- 
ceived a  great  deal  of  attention  since  it  first  confronted  Reynolds  [23] 
and  Taylor  [29]  . 

The  present  paper  opens  with  an  examination  of  the  physical  bases  for 
the  mathematical  description  of  turbulence  by  probability  distributions, 
and  for  the  clos\ire  problem  which  arises  from  this  descr:' _jtion.   After  dis- 
cussing these  matters,  we  shall  take  up  analytical  approaches  to  the  clo- 
sure problem  which  have  been  proposed  by  Deissler  [7:.8]^  Hopf  [lO,ll], 
Kraichnan  [l5,l6],  Proudman  and  Re  id  [21,22],  and  Tatsimil  [27,28].   We 
shall  describe  some  mathematical  inter-relations  and  then  Illustrate  cer- 
tain consequences  of  the  several  approaches  by  applying  them  to  an  ideal- 
ized turbulent  convection  problem  for  which  they  all  yield  explicit  solutions, 

2.    THE  PHYSICAL  NATURE  OF  TURBULENCE 

Fluid  flows  which  are  called  turbulent  by  fluid  dynamicists  typically 
have  two  related  characteristics:   a  complicated  variation  of  velocity  with 
position  and  time,  and  an  instability  to  small  disturbances.   In  experiments. 


many  kinds  of  small  disturbances  are  present  and  -uncontrollable.   Thus, 
the  instability  makes  it  impossible  to  predict  the  velocity  at  any  given 
time  and  place  from  a  knowledge  of  the  gross  initial  and  boundary  condi- 
tions. 

Let  us  consider  the  classic  example  of  water  which  flows  from  a 
reservoir  through  a  long  circular  pipe  into  a  second  reservoir  at  lower 
pressure.    Let  this  experiment  be  repeated  many  times  in  as  closely 
identical  fashion  as  possible.   If  the  flow  is  slow  enough,  it  is  fully 
laminar:   The  velocity  varies  smoothly  and  monotonically  from  wall  to 

axis  of  the  pipe,  and  it  is  very  nearly  reproducible  from  one  trial 

2 
(realization)  of  the  experiment  to  another.    If  the  flow  is  fast  enough, 

however,  the  velocity  in  a  typical  realization  displays  complicated  small- 
scale  variations  in  space  and  time,  and  these  variations  are  found  to 
change  unpredictably  from  one  realization  to  the  next.   At  a  representa- 
tive point  in  the  pipe,  the  measured  velocity  in  a  typical  realization 
will  differ  by  several  percent  from  the  mean  velocity  at  that  point  ob- 
tained by  averaging  over  an  ensemble  of  very  many  realizations.   More- 
over, the  velocity  fluctuations  remain  at  this  level  no  matter  how  hard 
we  try  to  keep  the  experimental  conditions  identical  to  all  realizations. 
We  shall  now  define  the  turbulent  velocity  field  as  the  difference  between 
the  actual  velocity  field  in  a  given  realization  and  the  ensemble-mean 
velocity  field.   Its  essential  characteristic  is  impredictability  under 


For  a  detailed  description  of  turbulent  pipe  flow,  see  Laufer  [JLTJ  and 
Townsend  [52]  . 

T!he   actual  experimental  difficulties  in  setting  up  accurately  reproduc- 
ible laminar  flows  can  be  very  severe. 


attainable  control  of  experimental  conditions. 

As  we  have  stated^  the  reason  for  the  appearance  of  the  tvurbulent 
velocity  field  appears  to  be  instability  to  small  disturbances.   These 
include  fluctuations  in  temperature  and  density  of  the  fluid,  velocity 
fluctuations  in  the  reservoir  (where  supposedly  the  fluid  is  at  rest), 
and  coupling  of  the  flow  to  vibrations  of  the  walls  of  the  pipe.   The 
dimensionless  parameter  which  principally  determines  whether  the  flow 
is  stable  or  not  is  the  Reynolds  number  DU/v,  where  D  is  the  pipe  dia- 
meter, U  the  mean  velocity  on  axis,  and  v  the  kinematic  viscosity.   For 
low  enough  Reynolds  number,  the  velocity  perturbations  in  the  pipe,  in- 
duced by  the  sources  just  cited,  are  damped  out  by  viscosity.   For  high 
enough  Reynolds  number,  the  viscous  friction  is  overbalanced  by  the  shear- 
ing action  of  the  existing  flow  upon  the  perturbation  velocity  [l8j  . 
This  action  causes  the  perturbation  velocity  to  draw  energy  from  the 
existing  flow  and  grow  rapidly  into  the  macroscopic  fluctuations  which 
characterize  the  turbulent  state.   The  process  is  demonstrated  dally  for 
those  of  us  who  stir  milk  into  our  coffee. 

The  pipe  flow  experiment  is  representative  of  actual  turbulent  shear 
flows  in  that  the  tiirbulent  velocity  at  any  point  is  typically  only  a  few 
percent  of  the  mean  velocity,  no   matter  how  high  the  Reynolds  number  is. 
Its  presence,  however,  is  associated  with  a  profound  change  in  the  charac- 
ter of  the  mean  flow.   In  the  fully  laminar  state,  the  mean  velocity  pro- 
file across  the  pipe  is  parabolic.   In  the  turbulent  state  at  high  Reynolds 
numbers,  the  mean  velocity  is  almost  constant  over  most  of  the  pipe  cross- 
section  and  falls  very  rapidly  to  zero  in  a  narrow  botindary  layer  at  the 
wall.   The  flow  resistance  is  very  high  compared  to  that  in  the  laminar 
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state^  and,  except  in  the  iimnediate  vicinity  of  the  wall;,  the  radial  trans- 
port of  momentum  is  dominated  by  the  turbulent  velocity  rather  than  by 
viscosity. 

The  unpredictability  and  complicated  form  of  the  turbulent  velocity 
field  lead  one  to  describe  it  by  probability  distributions  over  a  large  en- 
semble of  realizations.  Even   If  the  exact,  detailed  structure  of  a  single 
realization  could  be  measured,  it  vould  never  be  reproduced  in  any  repeti- 
tion of  the  experiment. 

A  paradoxical  feature  of  turbulence,  but  one  which  is  very  fortunate 
for  the  theorist,  is  that  extreme  sensitivity  of  individual  realizations 
to  small  perturbations  is  accompanied  by  an  insensitivity  of  the  probabil- 
ity distribution  to  the  nature  and  probability  distribution  of  the  pertur- 
bations.  Suppose  that  we  place  in  the  entrance  to  the  pipe  a  grid  of 
wires  which  substantially  disturbs  the  flow  and  thus  dominates  the  pertur- 
bations which  lead  to  turbulent  breakdown.   In  a  region  downstream  of  the 
grid,  the  probability  distribution  of  the  turbulent  velocity  will  depend 
on  the  parameters  of  the  grid.   Further  downstream,  however,  the  proba- 
bility distribution  will  be  almost  independent  of  the  nature  of  the  grid, 
and  of  its  presence,  provided  the  Reynolds  number  is  high  enough.  We 
shall  call  the  region  of  sensitivity  to  the  nature  of  the  triggering  dis- 
turbances the  transition  region,  and  the  region  of  insensitivity  that  of 
fully-developed  turbulence.   In  the  latter  region,  the  probability  distri- 
bution of  the  turbulent  velocity  is,  in  addition,  highly  insensitive  to 
small  disturbances  of  a  local  nature.   The  insensitivity  is  only  statisti- 
cal, however.   In  any  realization,  the  detailed  velocity  field  is  highly 
unstable  to  small  disturbances.   The  qualification  should  also  be  raised 


that  the  observed  statistical  insensltivity  might  be  less  striking  if 
more  detailed  and  accurate  statistical  measurements  could  be  made. 

5.   MATHEMATICAL  FORMULATION 

Let  us  write  the  total  velocity  field  In  a  realization  as 

(3.1)     v.(x,t)  =  v^(x,t)  +  u^(x,t), 

where  v  (x,t)  is  the  ensemble -average  of  v  (x,t)  and^  by  our  previous 
definition,  u  (x^t)  is  then  the  turbulent  velocity.   The  usual  statis- 
tical  description  of  u.(x,t)  is  by  the  infinite  set  of  moments 

\J^.{x,t;x',t')   =  <^u^(x,t)u^(x',t')^  , 


where  (     )   denotes  ensemble  average.   A  more  compact  description,  but 
one  which  so  far  has  not  been  much  used,  is  by  the  distribution  functions 
P-,(u,x,t),  P  (u,x,t;u' ,x',t' ),  etc.  defined  as  follows:  P.  (u,x,t)d  u  is 
the  probability  that  u(x,t)  lies  in  the  volume  d  u  about  the  point  u  in 
velocity  space,  P  (u,x,t;u' ,x', t ' )d  ud  u'  Is  the  joint  probability  that 

u(x,t)  lies  within  the  volume  d  u  at  u  and  u(x',t')  lies  within  the 

5 
voliime  d  u'  at  u',  etc.   Each  P  is  equivalent  to  an  infinite  subset  of 

the  moments  (5-2).   (Experimental  measurements  have  been  made  of  both 

moments  and  distribution  functions.) 

We  may  now  state  the  theoretical  problem  of  turbulence  as  follows: 

We  are  given  the  initial  mean  flow  v  (x, t  ),  together  with  the  complete 

set  of  moments  which  specify  the  distribution  of  the  initial  turbulent 


velocity  field |tu(x,t)   (if  the  latter  is  nonzero).   We  are  also  given 
a  complete  statistical  description_,  in  terms  of  moments,  of  whatever 
external  forces  or  perturbations  may  he  acting  during  the  course  of  the 
flow.   We  wish  to  find  v.(x,t);,  together  with  the  complete  set  of  moments 
(3. 2)  at  all  times  after  t  . 

Let  us  assume  that  the  fluid  is  incompressible,  of  unit  density,  and 
(apart  from  perturbation  effects)  obeys  the  Navler-Stokes  equation 

(5-3)       Q-t-^^\--\  ^-i:  ' 

V.        ■  mi 

where  v  is  a  constant  kinematic  viscosity  and  p  is  pressure.   The  pres- 
sure tenn  in  (5-5)  may  be  expressed  as  a  2nd -degree  polynomial  function- 
al of  the  velocity  field  by  use  of  the  boundary  conditions  and  the  diver- 
gence condition 

(5.^)     Sv^/Sx^  =  0. 

This  reduction  is  simplest  when  the  fl\ild  is  infinite  and  p  vanishes  at 
Infinity.   Then 

,    P        n   Sv  (x',t)   av  (x',t) 
(3.5)     P(x.t)=^     Ir^'l    -^ ^ ^^'' 


m         i 


and  (5.5)  may  be  written  in  the  form 


(5.6)      (n-^^O^i^-i^imn^^^^VJ^ 


where 


dx  Sx  /  dx    V   "        Sx.  Sx  Jbx 
1  my   n    V  i  n/   m 


and  V"^f(x)  =  -(4it)"^   /  |x-x' |"^f(x')d\'  for  any  f. 

We  have  so  far  not  represented  the  small  disturbances  which  play- 
so  important  a  role  in  turbulence.   In  order  to  include  their  effects 
accurately^  we  would  have  to  examine  in  detail  the  physical  mechanisms 
responsible.   Instead^  we  shall  be  guided  by  the  empirical  finding  that^ 
in  many  flows,  there  exists  a  fully-turbulent  state  which  is  insensitive 
to  the  precise  nature  of  the  disturbances.   The  simplest  assxamption  is 
to  lump  the  effects  of  all  disturbances  into  distortions  of  the  initial 
turbulent  velocity  field  and  then  assume  that  {3-3),   without  additional 
terms  and  without  perturbations  in  the  boundary  conditions,  holds  exactly 
throughout  the  course  of  the  flow.   This  implies  that  disturbances  during 
the  course  of  the  flow  merely  nudge  the  system  from  one  solution  of  (3'3) 
to  a  solution  corresponding  to  another  Initial  condition.   We  shall  adopt 
this  specialization  of  the  turbulence  problem.  We  must  recognize,  how- 
ever, that  it  cannot  be  expected  to  yield  an  accxirate  description  of  the 
transition  stage  of  turbulence.  Moreover,  its  adequacy  in  representing 
any  turbulent  state  is  only  a  conjecture  which  mxist   be  tested  by  its  con- 
sequences. 

In  principle,  we  may  now  solve  the  turbulence  problem  as  follows. 
We  assume  an  initial  mean  field  v(x, t  )  and  initial  values  for  the  set  of 
moments  (3-2) .   Then  we  construct  an  initial  ensemble  of  fields  v(x,t  )  + 


u(x, t  )  which  yield  these  moments.  We  solve  (5«5)  for  each  of  these 
initial  conditions  (assuming  solutions  exist  for  all  t)  and,  from  the 
solutions,  construct  v(x,t)  and  those  moments  (5«2)  which  are  of  in- 
terest. 

This  procedure  is  not  a  practical  one.   Few  solutions  of  (3.5)  are 
known  even  for  simple  initial  conditions.   To  carry  out  the  above  pro- 
gram accurately  for  a  large  set  of  complicated  initial  conditions  would 
Involve,  at  present,  a  prohibitively  large  numerical  computation  effort. 
Even  if  this  were  done,  there  would  remain  the  task  of  abstracting  from 
the  mass  of  calculations  an  insight  into  the  essential  dynamical  nature 
of  the  turbulence.   We  are  led.  Instead,  to  seek  a  method  of  solution 
which  takes  us  directly  from  the  Initial  set  of  moments  to  the  values, 
at  later  times,  of  the  moments  of  Interest.   To  do  this,  we  construct, 
from  (5.2),  equations  of  motion  for  the  moments  themselves.   Therein  lies 
the  closure  problem  which  is  to  be  our  principal  concern. 

h.        THE  CLOSURE  PROBLEM 

Let  us  now  specialize  to  an  infinite  fluid  for  which  (5.6)  is  valid. 
This  will  simplify  the  mathematics,  and  we  may  hope  that  the  essential 
nature  of  the  results  will  be  unchanged  in  the  more  general  case.   By 
averaging  (5.6)  over  the  ensemble,  we  obtain 


=  -  l^imn^^^Vil^^^^^^^)- 


V  (x,t)v  (x,t) 


Subtracting  (^.1)  from  (5.6),  we  obtain 


(^.2)      (^-  vv2)uJx,t).P.^(V)[^ 


Vjjj(x,t)u^(x,t) 


i  ^l.m(" 


u  (x,t)u  (x.t)  -  U  (x,t:x,t) 


where  we  note  P   (v)  =  P   (V) .   Now,  by  multiplying  (^.2)  with  u.(xSt), 

u.(x',t)u  (x",t"),  etc.  and  averaging,  we  obtain  the  infinite  set  of  equa- 
J  -~     r  ^^ 

tions 


^ 


^    V  V  ^)  U.  .(x,t;x',t')  +  P,   (V  ) 
X  y  ij  ~   -        imn^  X 


V  (x,t)U  .(x,t;x',t') 


(^.5) 


=  -  \  ^In^^V 


U   .(x,t;x,t;x',t') 
mnj  ~   ~   -w 


^-  ^^x7"ijr(^^>i^'^*'^i^"^*") 


+  P,   (V  ) 
imn  X 


V  (x,t)U  .  (x,t;x',t';x",t") 


-  P   fv  ) 
2  imn^  X'' 


U   .  (x,t:x,t:x',  t' 2x",t" ) 


U   (x,t;x,t)U.  (x',t';x",t") 
mn     ~.    jr  ~     -^ 


These  equations  are  formidable  in  appearance  (and  not  deceptively 
sol).   To  show  their  essential  structure,  we  shall  rewrite  (^.1)  and  (^.3) 
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in  the  schematic  form 


and 


(h.y) 


(Jt-   vv2)v+|p(V)[vv]   =   -i 


It  ■  "  "^  /  U^  +  P(V) 


^  -  V  V^  j  U^  +  P(V) 


V  U, 


V  U, 


2^2 


-  2  P(^)U3  ' 


-  \   P(V) 


\   -  U2U2 


where  V>     denotes  a  second-order  moment  and  so  forth.   We  see  that  (^.l) 
and  (^.5)  are  a  coupled  set  of  equations  that  link  moments  of  all  orders. 
In  general,  the  equation  for  the  nth-order  moments  involves  linearly  the 
(n  +  l)th-order  moments.   If  we  consider  the  integration  of  these  equa- 
tions forward  in  time,  it  is  clear  that  the  value  of  any  moment  at  later 
times  is  a  functional  of  the  initial  values  of  the  full  infinite  set  of 
moments. 

The  coupling  of  moments  of  all  orders  arises  from  the  right  side  of 
(4.2)  and  is  a  direct  consequence  of  the  nonlinearity  of  the  Navler-Stokes 
equation.   The  nonlinearity.  In  turn,  is  the  mathematical  embodiment  of 
the  self-stirring  of  the  velocity  field.  We  should  not  be  surprised  that 
the  moments  are  coupled.   In  any  realization  the  velocity  at  a  given  time 
and  place  is  determined  by  the  entire  initial  velocity  field.   It  is 
plausible,  then,  for  the  mean  velocity  at  that  place  to  depend  on  the  full 
initial  probability  distribution. 
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Our  previous  discussion  suggests,  however,  that  the  fiill  depen- 
dence of  given  moments  on  the  structure  of  the  initial  distribution 
is  far  more  information  than  we  want.   Indeed,  the  Idealization  of  the 
turbulence  problem  as  one  with  stochastic  initial  conditions  and  fully 
determinate  equations  of  motion  was  based  on  the  hope  that,  within  limits, 
the  evolution  of  the  turbulence  was  substantially  independent  of  the  de- 
tails of  the  initial  distribution. 

Let  us  therefore  specialize  our  problem  to  the  case  where  the 

initial  distribution  is  of  the  simplest  possible  form:  a  multivariate 

normal  distribution,  which  is  completely  specified  by  v  (x, t  )  and 

U, .fx.t  :x',t  ).   Eventually  we  would  like  to  be  able  to  solve  the  sta- 
ij  ~  0  —   0 

tistical  problem  for  a  variety  of  types  of  initial  distribution,  so  as 
to  see  just  what  sensitivity  to  this  choice  the  solutions  actually  dis- 
play.  The  multivariate  normal  distribution  seems  a  logical  starting 
point  because  it  is  simplest  and  beca\ise,  in  any  event,  we  do  not  really 
know  what  distribution  best  represents  the  physical  perturbations  which 
occur  in  experiments. 

At  this  point  we  come  to  a  fundamental  fact  which  is  distressing. 
Given  a  multivariate  normal  initial  distribution,  it  seems  clear  that  any 
given  moment  at  later  times  can  be  expressed  as  a  functional  of  v  (x, t  ) 
and  U  .(x,t  ;x' ,t    ).      If,  however,  we  attempt  to  find  the  explicit  form 


o 


f,  say,  v.(x,t)  by  integration  of  the  moment  equations  (^.l)  and  (^.5) 


1  '-- 


in  time,  we  obtain  an  infinitely  complicated  functional  of  v^(x,t  )  and 

U  fx,t  :x',t  ).   For  example,  if  we  do  the  integration  in  the  most 
ij  ~  0  ~   0 

straightforward  way,  by  expansion  in  a  Taylor  series  in  t-t^,  we  find 
that  V  (x,t)  is  an  infinite  functional  polynomial  series  in  these  quan- 
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titles.   Here  lies  the  essence  of  the  closure  problem.   We  siispect  on 
intuitional  grounds  that  a  simple  initial  distribution  should  result  in 
a  reasonably  simple  distribution  at  later  times.   Yet  we  find  that  the 
mathematical  relation  between  Initial  and  later  distributions  is  infinite- 
ly complicated. 

Perhaps  the  first  conclusion  suggested  by  this  situation  is  that  our 
formulation  is  inappropriate.   Perhaps  there  is  a  description^  other  than 
by  moments,  in  terms  of  which  simple  initial  distributions  would  give 
manifestly  simple  later  distributions.   Such  a  formulation  has  not  been 
reported^  and  there  is  an  elementary  reason  why  one  should  be  hard  to 
find.   A  simple  initial  distribution,  such  as  the  normal  one  we  have  adopt- 
ed_,  does  not  mean  a  simple  initial  field  u  (x,t  )  in  typical  realizations. 

Even  if  we  make  the  further  requirement  that  U^  .(x, t^tx' .t^)  be  a  smooth 

^  ij  —  0  --   0 

and  simple  function,  a  typical  realization  of  u  (x,t  )  will  be  a  compli- 
cated  function  of  x.   The  solution  of  the  Navier-Stokes  equation  for  such 
realizations  is  very  difficult.  And,  since  the  various  realizations  do 
not  interact,  the  task  is  not  made  easier  by  considering  an  ensemble  in- 
stead of  a  single  system.  As  long  as  our  formulation  of  the  statistical 
problem  is  based  ultimately  on  the  Navier-Stokes  equation,  the  quest  for 
a  simple,  exact  statistical  solution  would  seem  to  depend  on  finding  a 
simple  general  solution  of  the  Navier-Stokes  equation  itself.  Let  us 
hope  that  there  is  a  flaw  in  this  argument i 

At  this  juncture,  it  is  illiuninating  to  compare  the  turbulence  prob- 
lem with  that  of  the  Gibbsian  statistical  mechanics  of  Hamiltonian  sys- 
tems [^51]  .   The  Gibbs  method  yields  an  explicit  distribution  function 
for  a  Hamiltonian  system  in  thermal  equilibrixim.  To  obtain  this  result. 
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one  need  not  solve  the  dynamical  equations.  We  may  ask  vhy  a  similar 
program  cannot  be  carried  out  for  the  turbulence  problem. 

The  answer  appears  to  be  tvo-foldi  First,  the  systems  which  can  be 
treated  successfully  by  the  Gibbs  method  are  special  ones.   Second,  even 
for  these  systems,  the  information  obtained  in  closed  form  is  actually 
a  small  part  of  a  complete  statistical  description;  it  does  not  include 
some  types  of  information  which  are  of  interest  in  turbulence  theory. 
The  special  property  which  is  central  to  the  Gibbs  theory  is  the  exist- 
ence of  Liouville's  theorem.   The  latter  states  that  the  motion  of  the 
"phase  cloud"  representing  an  ensemble  is  incompressible  in  the  phase 
space  of  the  system;  that  is  to  say,  it  preserves  measure.   This  pro- 
perty depends  on  the  conservative  character  of  the  system.   In  turbu- 
lence theory,  the  phase  space  is  the  function  space  of  v,  a  correspon- 
dence which  has  been  elegantly  formulated  by  Hopf  [loj .  However,  dis- 
sipation is  important  in  turbulence,  and,  as  a  consequence,  measure  in 
the  phase  space  is  not  even  approximately  conserved.   In  the  presence 
of  dissipation,  the  Gibbs  method  does  not  lead  to  closed  expressions 
for  distributions,  and  it  remains  necessary  to  solve  the  equations  of 
motion. 

The  second  part  of  the  answer  is  that,  even  for  conservative  sys- 
tems in  thermal  equilibrium,  the  Gibbs  method  gives  only  the  simultane- 
o\is  distribution  in  closed  form.   To  obtain  moments  with  non-simultane- 
ous time-arguments,  one  again  must  return  to  the  equations  of  motion. 
If  the  latter  are  nonlinear,  there  arises  a  closure  problem  quite 
analogous  to  that  of  turbulence  theory.   Subordination  of  the  role 
accorded  non- simultaneous  averages  is  not  very  appropriate  in  turbulence 
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theory  because  such  averages  give  valuable  information  about  the  dynamics 
of  the  nonlinear  interaction.  Moreover^  transport  phenomena,  which  always 
occur  when  there  is  dissipation,  are  most  naturally  expressed  by  non-simul- 
taneous averages:  Transport  implies  that  something  goes  from  a  place  at  a 
given  time  to  another  place  at  a  later  time . 

5.   SCBEMES  FOR  APPROXIMATE  CLOSURE 

The  unavailability  of  a  simple  method  of  following  exactly  the  evolu- 
tion of  a  simple  initial  velocity  distribution  suggests  that  approximations 
be  sought.   This  seems  particularly  appropriate  because  approximate  empiri- 
cal observations  are  the  real  basis  for  the  intuition  that  simple  initial 
distributions  ought  to  lead  to  reasonably  simple  later  distributions. 

There  have  been  two  broad  classes  of  attempts  to  close  off,  or  circum- 
vent, the  infinite  set  of  moment  equations  by  means  of  approximations.   The 
first  consists  of  essentially  phenomenological  approaches  in  which  statisti- 
cal properties,  or  models  of  dynamical  processes,  are  postulated  on  physical 
grounds.   Examples  are  the  universal  equilibrium  theory  of  Kolmogorov  [^,12, 
15],  Heisenberg's  eddy-viscosity  theory  [^^9]^  and.  the  maximum-dissipation 
theory  of  Malkus  [19] .   The  second  class  consists  of  schemes  for  systematic 
analytical  approximation.   Our  basic  interest  in  the  present  paper  will  be 
in  the  mathematical  structure  of  turbulence;  consequently,  we  shall  not 
discuss  the  first  class  of  approaches.   For  proper  perspective,  however,  it 
is  highly  important  to  note  that  almost  all  progress  to  date  in  the  physical 
understanding  of  turbulence  has  come  from  phenomenological,  heuristic  theo- 
ries such  as  the  three  cited.   It  is  hardly  unfair  to  say  that,  so  far,  the 
principal  achievement  of  the  analytical  approaches  has  been  to  clarify  as- 
pects of  their  own  structures  and  shortcomings. 
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5-1-   Isotropic  turbulence 

The  analytical  approximations  have  been  carried  out  principally 
for  the  special  case  of  isotropic  or  homogeneous  turbulence.  We  shall 
treat  them  in  this  context,  but  we  want  first  to  discviss  the  degree  of 
relevance  this  artificial  problem  has  to  the  understanding  of  actual  tur- 
bulent flows.   The  internal  dynamical  processes  which  act  on  the  turbulent 
velocity  field  u(x,t)  are  three:  viscous  friction,   interaction  of  mean 
field  with  turbulent  field,  and  self-interaction  of  the  turbulent  field. 
They  are  represented  in  (^.2)  by  the  terms  involving  v,    \y   (x,t)u  (x,t)J, 
and  [u  (x^,t)u  (x,t)  -  U  (x,t;x,t)],  respectively.   Only  the  last  of  the 
three  processes  gives  rise  to  the  closure  problem,  as  we  have  already 
noted.   It  is  therefore  natural  to  try  to  understand  the  problem  in  the 
context  of  isotropic  turbulence:   There  is  no  mean  field,  and  the  self- 
interaction  of  the  turbulent  field  provides  the  only  dynamical  process 
besides  the  essentially  simple  one  of  viscous  friction.  Moreover,  geo- 
metrical complications  are  minimum  for  this  case. 

There  are  two  dangers,  however,  in  concentrating  on  the  isotropic 
case.   The  first  is  one  of  perspective:   Even  if  the  closure  problem 
were  solved  to  complete  satisfaction,  there  would  still  remain,  for  shear 
turbulence,  the  very  formidable  task  of  treating  the  interaction  of  mean 
field  and  turbulent  field.   The  magnitude  of  this  task  may  be  illustrated 
by  the  observation  that  its  most  elementary  form  (valid  only  for  Infini- 
tesimal turbulence  levels)  comprises  the  laminar  stability  theory  of  shear 
flow  [18] .   The  second  danger  is  that,  despite  the  geometrical  simplicity, 
the  problem  of  finding  adequate  closure  approximations  may  actually  be 
more  difficult  for  isotropic  turbulence  than  for  shear  flows.   In  a  shear 
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flow,  important  aspects  of  the  overall  dynamics  may  be  dominated  "by  the 
mean  field-turbulent  field  interaction.   If  so^  cruder  approximations 
to  the  turbulent  field  self-interaction  may  be  permissible  than  in  the 
Isotropic  problem,  where  this  self -interaction  is  all-important. 

With  this  preface,  let  us  specialize,  in  the  remainder  of  this  Sec- 
tion, to  homogeneous,  isotropic  turbulence,  for  which  v  (x,t)  =  0.  Then 
(4.2)  reduces  to  the  form 


(5-1)      Qt   -   ^   ^')  \^^'^^    -   -  -2   ^imn(^) 


u  (x,t)u  (x,t) 


Let  us  assume,  as  before,  that  u.(x,t  )  is  multi-variate  normal.   Let  v 
be  the  initial  root-mean-square  value  of  any  velocity  component,  and  let 
I     be  a  characteristic  correlation  length  (macro-scale  [hj  )    of  the  ini- 
tial covariance  U.  .(x,t  :x',t  ).   In  order  to  exhibit  clearly  the  rela- 
tive magnitudes  of  the  terms  in  (5-1),  we  shall  measure  lengths  and 

p 
velocities  in  units  I      and  v  ,  and  times  in  the  unit  I      /v,  which  is  the 

characteristic  time  for  viscous  decay  of  a  velocity  structure  of  scale  I    . 

Using  the  new  units,  which  we  shall  adopt  for  the  remainder  of  Section  5, 

we  may  rewrite  (5-1)  in  the  form 

(5-2)     Q^  -    V^)  u^(x,t)  =  -  I  RoPi^(V)  [u,U.t)u^(x,t)]  , 
where 

(5-5)    Ro  =  Vo/^ 

is  the  characteristic  Reynolds  number  for  the  initial  state. 
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^.2  Expansion  in  powers  of  Reynolds  number 

A  formal  solution  of  (5.2)  may  be  obtained  by  writing 


(5.M     u^(x,t)  =  ^     R\^^''kx,t), 


n=0   0 


substituting  in  (5.2),  and  equating  coefficients  of  like  powers  of  R  on 
each  side  of  the  equation.   In  this  way  we  find 


(0).    .^  _  rk.^._.  ^1-5/2 


(5.5)     u/"^(x,t)  =  Un(t-t 


1 


0' 


exp 


-lx-x'|^A(t-tQ)lu^(x',tQ)d\s 


(n), 


and,  by  iteration,  we  obtain  u.    (x,t)  as  a  functional  power  of  u  (x,t  ) 
of  degree  n  +  lo   The  expansion  (5.*^-)  then  permits  us  to  write  any  moment 

of  u  (x,t)  as  an  infinite  series  of  linear  functionals  of  the  moments 

1  - 

U^.(x,t  :x',t  ),  U^.Jix,t^;x' ,t^;x" ,t^),    etc.    These  moments,  in  turn, 
all  may  be  expressed  in  terms  of  U  .(x,t  :x',t  )  according  to  the  assump- 
tion that  the  initial  distribution  is  normal.   In  this  way,  we  obtain  a 
complete  formal  solution  of  the  statistical  problem.   It  illustrates  our 
earlier  statement  that  the  exact  expression  of  moments  at  later  times  in 
terms  of  initial  moments  is  one  of  infinite  complexity. 

The  most  straightforward  sequence  of  closure  approximations  for  the 
present  problem  consists  in  truncating  the  iteration  expansions  for  the 
moments  at  successively  higher  powers  of  R  .   The  fundamental  limitation 
to  this  scheme  is  that  we  may  reasonably  expect  the  approximations  to  con- 
verge rapidly  only  if  R  «  1,  in  which  case  the  nonlinear  term  represents 
a  small  perturbation  to  (5.2)  in  a  typical  realization.   This  condition  is 
satisfied  during  the  final  period  of  decay  of  turbulence  produced  behind 
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a  grid.   For  this  case,  the  zeroth  term  in  the  Reynolds  number  expansion 
(i.e.,  pure  viscous  decay)  has  been  studied  by  Reissner  and  by  Batchelor 
and  Townsend  [h~\  ,    who  obtained  interesting  results  concerning  the  asymp- 
totic distribution.   Several  higher  terms  in  the  expansion  for  U. .  and 
U  .   have  been  investigated  by  Deissler  [7^8]  . 

For  R   »  1,  which  is  a  ca.se  of  great  interest,  we  cannot  expect 
the  iteration  expansion  to  yield  adequate  approximations.   There  is  an 
exception,  however.   For  t-t   sufficiently  small  compared  to  character- 
istic times  of  the  turbulence,  it  is  reasonable  that  the  expansion  should 
converge  well  even  for  R  »  1.   To  see  this,  we  note  from  (5.2)  that 

\h   u.(x,t)/dt  J     must  be  a  polynomial  in  R  of  degree  n.   Thus,  only 
1  -        t=tQ  0 

low  powers  of  R  appear  in  low  orders  of  a  Taylor  expansion  in  t-t  .   A 

similar  situation  will  be  examined  in  Section  6.2. 

Although  the  iteration  expansion  can  be  expected  to  yield  useful 
approximations  only  for  a,  quite  limited  range  of  parameters,  it  is  of 
great  value  as  an  exact  formal  solution  against  which  to  compare  approxi- 
mate solutions  obtained  by  other  means.   As  we  shall  point  out  shortly, 
several  other  systematic  approximation  sequences  which  have  been  proposed 
may  be  represented  formally  as  summations  of  contributions  selected, 
according  to  well-defined  rules,  from  all  orders  of  the  iteration  expan- 
sion. 

An  interesting  question  at  this  point  is  whether  an  expansion  in 
powers  of  1/R^  might  be  useful  for  R  »  1.   It  seems  doubtful  that  such 
a  procedure  would  work.   The  reason  is  a  fundamental  asymmetry  in  the 
roles  of  the  nonlinear  forces  and  viscous  forces.   For  R  «  1,  the  non- 
linear forces  are  a  small  perturbation.   But  for  R  »  1,  the  viscous 
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forces  are  not  a  small  perturbation:   The  half-life  for  viscous  decay  of 
the  turbulent  kinetic  energy  is  the  order  of  the  characteristic  time  ^^/v 
(original  units)  for  nonlinear  effects  [U] .   This  behavior  is  intimately 
related  to  the  fact  that  solutions  of  the  inviscid  Navier-Stokes  equation 
(R  =  oo )  can  develop  singular  shear-fronts  in  a  finite  time,  even  if  the 
initial  conditions  are  healthy. 


^♦3   Expansion  in  cumulants 

As  we  have  seen,  the  iteration  expansions  for  the  moments  may  be 
derived  from  a  similar  expansion  for  the  velocity  field  in  each  individual 
realization.   For  this  reason,  one  suspects  that  the  iteration  approxima- 
tions do  not  exploit  fully  the  fact  that  we  are  interested  only  in  averages. 
The  approximations  we  shall  consider  next  differ  in  this  respect.   They  are 
intrinsically  statistical,  and  they  have  no  meaning  as  approximations  on 
the  individual  realizations. 

In  19^1,  Millionstchikov  [20]  introduced  the  assumption  that  fourth- 
order  moments  of  the  distribution  of  simultaneous  velocity  amplitudes  are 
related  to  second-order  moments  as  in  a  normal  distribution.   This  quasi- 
normality  approximation  has  been  applied  to  isotropic  tiirbulence  dynamics 
by  Proudman  and  Reid  [21]  and  Tatsumi  [27,28].   A  similar  approximation  for 
the  two-time  velocity  distribution  has  been  exploited  by  Chandrasekhar  [6] 
and  Heisenberg  [9] .   Hopf  [ll]  has  pointed  out  that  the  quasi-normality 
approximation  is  equivalent  to  retaining  only  the  initial  terms  in  a  Gram- 
Charlier  expansion  for  the  velocity  distribution,  and  he  has  suggested 
that  a  systematic  approximation  sequence  could  be  based  on  this  expansion. 
We  shall  discuss  here  an  approximation  sequence  based  on  expansion  of  the 
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velocity  moments  into  ciimulants.   This  is  closely  related  to  the  Gram- 
Charlier  expansion. 

Let  P(y)  be  the  probability  distribution  for  a  single  random  vari- 
able y  which  has  zero  mean.   Then 


(5.6) 


P(y)dy  =1,       (y)  =   r  yP(y)dy  =  0 


The  cumulants  of  the  distribution  may  be  defined  as  the  coefficients  c 

r 

in  the  expansion 


(5.7)     $(a)  =  exp  -  |(y  /a  +  2Li     c^(ia)VnJ 


n=5 


where 


(5.8)     ^(a)  =  (exp(iay))    =   /  exp(iay)P(y)dy 

-  00 

is  the  characteristic  function.   Upon  evaluating  the  moments  by  use  of 
the  relation 


(5.9)  (y'')     =    (-i)^d^$(a)/d..^_ 


a^O 


we   find 


(5.10) 


-y 


=   3</)^ 


=   10 


</> 


+  C),^ 


C5   +    c^. 


15  (y^)^  +  loc^^  + 


^5  </) 


■=1*  *  %• 
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The  significance  of  the  cumulants  is  that  they  measure  the  degree  and  com- 
plexity of  the  deviation  of  P(y),  and  its  moments,  from  a  normal  distribu- 
tion. 

If  an  expansion  analogous  to  (5"7)  is  made  for  the  characteristic 
functional  [lO,llj  of  the  one-time  velocity  distribution,  a  sequence  of 
approximations  to  the  various  moments  can  be  formed  by  truncating  the 
expansion  so  as  to  admit  successively  higher  cumulants.   In  the  first  non- 
trivial  approximation,  we  admit  third-order  cumulants  only.   Then  [in 
correspondence  to  the  second  of  relations  (5.I0),  with  c,  =  6]   we  find 

(5.11)   U^j^(x,t;x',t;x",t;x"',t)  =  U^^(x,  t;x' ,t)U^(x",  t;x '"  ,t) 

+  U^Jx,t;x",t)U^^(x',t;x"',t)  +  U.^(x,t;x '"  ,t)U^^(x' ,  t;x"  ,  t) , 

which  is  the  quasi -normality  approximation  employed  by  Proudman  and 
Reid  [21]  and  Tatsumi  [27,28].   More  generally,  by  cutting  off  the  expan- 
sion after  the  cumulant  of  nth  order,  we  obtain  a  closed  expression  for 
moments  of  order  n  +  1  in  terms  of  moments  of  lower  orders.  We  shall 
call  the  sequence  of  approximations  so  obtained  the  cumulant-discard 
approximations . 

Equation  (5.11)  may  be  used  to  close  the  moment  equations  as  follows. 
From  (5-5)  we  form  equations  for  the  total  time-derivatives  of 
U.  .(x,t;x',t)  and  U.   (x, t;x' ,t;x" ,t) .   We  then  eliminate  the  fourth-order 
moments  from  these  equations  by  (5. 11),  thereby  obtaining  a  closed  set  of 
sim\iltaneous  equations  for  U.  .(x,t;x',t)  and  U.  ..  (x,  t;x' ,t;  x"  ,  t) .   The 
Third-order  moments  can  be  eliminated  from  this  set  by  integration  so  that, 
finally,  we  end  up  with  an  Integro-differential  equation  for  U.  .(x,t;x',t) 
alone . 
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At   the   next   level,    ve   eliminate   fifth-order  moments  by  means   of  the 
second  cumulant-dlscard  approximation    ^corresponding  to  the   third  of  rela- 
tions   (5.10),    vith   c     =   0I .      We   thereby  obtain  simultaneous   equations   for 

5 

second-,  third-,  and  fourth-order  moments.   Then  we  eliminate  the  fourth- 
order  moments  by  integration  and  finally  obtain  closed,  simultaneous  equa- 
tions for  U.  .(x,t;x'  ,t )  and  U   (x, t;x' , t;x" , t) .   As  we  go  up  in  this  se- 
quence  of  approximations,  the  complexity  of  the  final  set  of  equations  which 
must  be  solved  increases  very  rapidly.   The  isotropy  condition  results  in 
substantial  simplifications,  however. 

The  cumulant- discard  approximations  bear  a  definite  relation  to  the 
iteration  expansion.   If  we  solve  the  final  equations  for  the  nth  cumulant- 
discard  approximation  by  expansion  in  powers  of  R  ,  we  find  that  the  expan- 
sion for  U   (x,t;x'  ,t)  contains  all  (even)  powers  of  R  .   It  agrees  pre- 
cisely  with  the  exact  iteration  expansion  to  order  R    .     This  is  a 
consequence  of  the  fact  that  the  exact  iteration  expansion  of  the  cumulant 
of  order  n  +  2  starts  with  H^ .      Thus,  the  neglect  of  this  cumulant  can 
affect  only  terms  above  a  certain  order  in  the  expression  for  U   .   For 
powers  of  R   greater  than  2n,    the  cumulant -discard  result  for  U   is  no 
longer  exact,  but  the  coefficient-function  of  each  power  of  R  bears  a 
definite  relation  to  that  in  the  exact  expansion.   The  coefficient  in  the 
exact  expansion  consists  of  a  sum  of  individual  terms,  and  that  in  the 
cumulant-discard  result  consists  of  only  some  of  these  terms,  selected 
according  to  well-defined  and  simple  rules.   We  shall  discuss  this  relation 
in  more  detail,  for  a  simpler  problem,  in  Section  6.5 • 

Since  the  nth  ciMiulant- discard  approximation  to  U. .  agrees  with  the 

J- J 

exact  iteration  expansion  to  order  R   , 
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it  is  clear  that  for  R^  «  1,  or  for  t-t  sufficiently  small,  this  approxi- 
mation will  be  indistinguishable  from  the  nth  iteration  approximation.   Of 
more  interest  is  the  behavior  of  the  cumulant-discard  approximations  for 
Rq  »  1.   The  fact  that  they  contain  all  powers  of  R  offers  the  hope  that 
they  may  converge  rapidly  when  an  expansion  in  R  does  not.   It  is  not  clear, 
at  the  present  time,  whether  this  is  the  case  or  not.  When  R  is  large,  it 
seems  unlikely  that  all  the  cumulants  can  be  small,  but  this  does  not  exclude 
the  possibility  that  their  magnitude,  or  dynamical  importance,  may  decrease 
with  increasing  order.   If  this  does  happen,  with  sufficient  rapidity,  it  is 
reasonable  to  expect  good  convergence.   The  existing  experimental  evidence 
on  the  magnitude  of  cumulants  [kj    is  difficult  to  interpret.   This  is  partly 
because  of  low  accuracy,  but  the  principal  reason  is  that  the  appropriate 
quantitative  theoretical  criteria  for  "smallness"  have  not  been  worked  out. 
Apart  from  the  question  of  the  quantitative  accuracy  of  the  cumulant- 
discard  approximations,  there  is  a  danger  that  they  may  lead  to  unphysical 
consequences  of  a  disconcerting  kind  at  high  R  .   This  arises  because  the 
distribution  functionals  obtained  by  truncation  of  the  expansion  in  cumu- 
lants are,  in  general,  not  positive-definite  in  the  function-space  of  u(x,t); 
that  is  to  say,  the  cumulant-discard  approximations  can  lead  to  negative 
probabilities.   This  seems  actually  to  occur  in  the  first  approximation. 
Proudman  and  Reid  [2l]  have  been  able  to  integrate  exactly  the  equation  for 
the  skewness  factor  S  =  \u   / /<fu   y  '    in  this  approximation.   Their 
result  shows  that  S  takes  values  which,  according  to  inequalities  derived 
by  Betchov  [5] ,  are  incompatible  with  a  positive-definite  distribution 
having  zero  fourth-order  cumulants. 
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The  generation  of  negative  probabilities  raises  the  possibility  that 
the  cumulant -discard  approximations  may  lead  to  non-positive-definite  wave- 
number  spectra  E(k)  for  the  kinetic  energy.   In  inhomogeneous  turbulence, 
they  may  lead  to  negative  kinetic  energy  density/— ru(x,t)l  yin  some  regions 
of  space.   The  cumulant-discard  approach  so  far  has  been  explicitly  vorked 
out  only  for  the  first  approximation  and  only  for  isotropic  turbulence. 
The  equations  in  this  case  have  been  integrated  to  give  the  evolution  of 
E(k)  only  for  very  short  times  or  small  R  ,  where  the  validity  of  the 
approximations  is  not  in  doubt.   Thus,  the  questions  we  have  raised  are 
unanswered  at  present.   However,  as  we  shall  show  in  Section  G.^,    the  cumu- 
lant-discard approximations  for  the  one-time  distribution  do  give  unphysical 
behavior  when  they  are  applied  to  the  convection  of  a  passive  scalar  field 
\|r(x,t)  by  a  random  velocity  field.   They  lead  there  to  a  situation  in  which 

^[jf(x,t)]   )  ,  which  is  analogous  to  mean  kinetic  energy  density,  goes  nega- 

5 

d  X,  which  should  be  conserved,  increases 


tive  in  places  and 


<[t(x.t)]2) 


without  limit.   The  difficulties  with  negative  probabilities  become  demon- 
stx-ably  severe  in  turbulence  dynamics  proper,  if  the  cumulant-discard  approach 
is  extended  to  the  many-time  distribution  so  as  to  give  non-simultaneous  as 
well  as  simultaneous  averages  [l4_]  .   The  existence  of  all  these  troubles  is 
discouraging  because  the  cumulant-discard  approach  seems  to  be  a  very  natural 
and  logical  one.   It  has,  moreover,  a  record  of  successes  in  other  physical 
problems .  . 

^.4.   Approximation  by  stochastic  models 

An  alternative  to  the  iteration  and  cumulant-discard  approaches  has 
recently  been  proposed  by  the  author  [l6j  .   The  basic  idea  is  to  replace  the 
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true  dynamical  problem  by  a  model  in  which  fictitious  couplings  are  intro- 
duced among  an  infinite  collection  of  flow  systems .  The  true  problem  then 
is  recovered  as  the  limit  in  which  all  the  couplings  vanish.  The  phaslngs 
of  the  couplings  change  at  random  as  we  go  from  one  individual  system  in 
the  collection  to  another.  Thus,  the  model  involves  a  new  stochastic  ele- 
ment which  is  absent  in  the  true  problem. 

The  additional  randomness  plays  a  paradoxical  role.   It  makes  the 
detailed  dynamics  of  the  collection  of  coupled  systems  infinitely  more 
complicated.   But  it  vastly  simplifies  the  moment  equations  obtained  by 
averaging  over  an  ensemble  of  realizations  of  the  collection.   This  is 
because  cancellation  effects,  arising  from  the  random  phasings  of  the 
couplings,  eliminate  many  of  the  complexities  which  the  nonlinearity  of 
the  dynamical  processes  would  otherwise  produce.  As  a  result,  the  model 
problem  leads,  without  any  approximation,  to  closed  equations  for  moments 
of  low  order,  In  contrast  to  the  infinite  set  of  coupled  equations  which 
characterize  the  true  problem.   We  note.d  in  Section  h-   that  the  considera- 
tion of  an  ensemble  could  not,  in  itself,  be  expected  to  alleviate  the 
fundamental  difficulty  of  solving  the  Navler-Stokes  equation  in  the  realiza- 
tions which  make  up  the  ensemble.   In  the  model  method,  we  circumvent  this 
limitation  by  making  many  systems  interact  with  each  other  in  a  stochastic 
fashion. 

Successive  approximations  in  the  stochastic  model  scheme  must  be 
sought  by  constructing  a  sequence  of  models  in  which  the  random  dynamical 
couplings  are  systematically  reduced.   This  must  be  done  in  such  a  fashion 
that  exact  closed  equations  (which  will  be  of  Increasing  complexity)  result 
at  each  stage.   So  far,  only  one  model  has  been  explicitly  constructed.  A 
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putative  second  model  has  been  investigated^  but  its  self- consistency 
has  not  been  proved.   At  present,  ve  cannot  be  sure  that  there  actually 
exists  an  infinite  sequence  of  models  vhich  converge  to  the  true  problem. 

The  stochastic  model  scheme  has  the  advantage  that  negative  probabil- 
ity troubles,  such  as  negative  spectra  or  energy-densities,  cannot  arise. 
This  is  because  the  final  equations  are  the  exact  description  of  the  dyna- 
mical model  and  do  not  involve  an^y  approximation  to  its  probability  distri- 
bution.  A  substantial  price  is  paid  for  this  advantage,  however.   The 
first  approximation  in  the  cumulant-discard  sequence  for  isotropic  turbu- 
lence leads  to  a  closed  equation  for  U.  .(x,t;x',t),  the  covariance  tensor 
for  equal  times.   In  the  first  stochastic  model  also,  closed  equations  are 
obtained  for  the  coveriance  tensor.   However,  it  is  necessary  to  determine 
the  full  non-simultaneous  tensor  U. .(x, t;x' , t ' ) .   The  simultaneous-argument 
problem  cannot  be  separated  off.   At  the  same  time  it  is  necessary  to  solve 
equations  which  determine  the  average  response  of  the  system  to  infinitesi- 
mal disturbances.   This  additional  information  is  actually  of  considerable 
interest  (cf.  Sections  2  and  ^),  but  the  equations  necessarily  are  more 
complex. 

If  the  result  for  U.  .(x,t:x',t)  in  the  first  stochastic  model  is 
expanded  in  powers  of  R  ,  a  definite  relation  is  found  to  the  exact  itera- 
tion expansion  for  U.  . (x,t;x' ,t) .   In  common  with  the  first  cumulant-discard 

ij  ~~   ~- 

approximation,  the  model  result  contains  all  (even)  powers  of  R  and  is 

2 
identical  with  the  iteration  expansion  to  order  R   .   The  coefficients  of 

the  higher  powers  of  R  represent  a  selection  of  terms  from  the  correspond- 
ing orders  of  the  iteration  expansion,  but  a  more  comprehensive  one  than 
that  included  in  the  first  cumulant-discard  approximation.   The  selection 
rules  are  given  in  [l6j . 
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When  applied  to  isotropic  turbulence,  the  first  stochastic  model  may 
he  given  an  alternative  interpretation  which  does  not  involve  the  coupling 
of  a  collection  of  flow  systems.   It  may  be  regarded  as  an  approximation 
on  the  dynamical  interaction  of  the  spatial  Fourier  amplitudes  of  the 
velocity  field  in  a  single  (infinite)  system  [15]  .   The  approximation  con- 
sists of  retaining  only  that  contribution,  to  the  phase  correlation  of 
each  triad  of  Fourier  amplitudes,  which  is  induced  by  the  direct  dynamical 
interaction  of  the  three  amplitudes. 

The  model  gives  a  description  of  isotropic  turbulence  dynamics  which 
appears  to  be  qualitatively  acceptable  [15J .   It  predicts  separate  energy- 
containing,  inertial,  and  dissipation  ranges  at  high  Reynolds  numbers.   In 
addition,  it  yields  plausible  statistical  stability  properties.   However, 
the  asymptotic  spectrum  law  for  the  inertial  range  disagrees  slightly  with 
that  of  the  Kolmogorov  theory,  and  it  is  probably  inaccurate  to  that  extent. 
It  is  not  yet  known  how  the  inertial  range  structure  is  altered  in  the 
tentative  second  model. 

6.   AN   APPLICATBDN  OF  THE  CLOSURE  APPROXIMATIONS 

6.1  Convection  by  a  random  velocity  field 

We  shall  now  illustrate  some  properties  of  the  several  closure  schemes 
by  applying  them  to  a  problem  for  which  they  all  yield  explicit  results 
that  can  be  compared  with  the  exact  solution.   For  this  purpose,  we  consi- 
der the  convection  of  marked  fluid  particles  by  a  prescribed  random  velocity 
field.   Since  the  initial  analytical  treatment  of  this  problem  by  Taylor  [29]^ 
it  has  received  extensive  study.   The  formulation  we  shall  employ  here  follows 
that  of  Batchelor  [l-j]  • 
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Consider  an  ensemble  of  incompressible  flows  in  which  the  velocity 
field  u(x)  has  an  isotropic,  multivariate  normal  distribution,  and  is 
independent  of  time  in  each  realization.   The  distribution  then  is  wholly 
characterized  by  the  covariance 


(6.1)  U,  .(x-x')  =  (u  (x)u  (x')) 


A   physical  embodiment  of  a  typical  realization  could  be  a  vessel  of  water 
agitated  at  low  Reynolds  number  by  a  very  large  number  of  small;,  steady 
stirring  devices^  located  and  oriented  at  random  throughout  the  fluid. 
(The  stirring  devices  might,  for  example,  be  ciliated  micro-organisms.) 
The  concentration  of  marked  particles  \)f(x,t)  obeys  the  equation 


(6-2)        U-^'>-"i^ 


where  K   is  the  molecular  diffusivity.   We  shall  treat  the  simplest  statis- 
tical problem  associated  with  this  equation:   Given  the  initial  mean  con- 
centration ^\|/(x,0))  ,    with  \|/(x,0)  statistically  independent  of  u,  find 
the  later  mean  concentration  <^ijf(x,t)^  .   Since  (6.2)  is  linear,  this 
problem  reduces  to  that  of  finding  the  average  Green's  function,  defined 

(6.5)  G(x-x',t)  =  (^(x,t))  ,       1r(x,0)  =  b^ix-x'), 

where   \lr(x,t)    obeys    (6.2)    and  5    (x-x' )    is   the   three-dimensional  Dirac   func- 
tion.      G{x-x' ,t-t')d  X  is   the  probability  that   a  marked  particle   introduced 
at   (x',t')    is    found  in  the  volume   d  x  at    (x,t). 
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If  we  multiply  (6.2)  by  various  powers  of  ^   and  average,  we  find 
a  closure  problem  quite  similar  to  that  of  turbulence  dynamics.   This 
occurs  because  (6.2)  is  nonlinear  in  stochastic  quantities,  although  it 
is  linear  in  the  dynamic  variable.  We  shall  apply  to  the  determination 
of  G(x, t)  all  three  of  the  closure  schemes  discussed  in  Section  5' 
First,  we  wish  to  establish  some  features  of  the  exact  solution  by 
appealing  to  simple  analytical  and  physical  arguments. 

As  in  Section  5,  let  I     and  v  be  the  correlation  length  and  root- 
mean-square  velocity  component  associated  with  u(x).   The  parameter 
which  plays  the  role  of  R  is 

(6.M  B^  =  ^^v^A.   5 

We  shall  confine  our  attention  to  the  case  B»  1.   (We  shall  not  adopt 
special  time,  length,  and  velocity  scales  here  as  we  did  in  Section  5  ) 

Two  asymptotic  ranges  may  be  distinguished  in  which  G(x^t)  has  a 
simple  form.   They  are  t  «  ^  /v  and  t  »  ^q/^q-      In  the  first  range^ 
the  distance  from  the  space  origin  travelled  by  almost  all  particles  is 
«  I    .   Thus,  they  have  very  nearly  their  initial  velocity  u(0).   The 
distribution  is  then  the  resultant  of  two  independent  Gaussian  processes; 
molecular  diffusion  and  a  Gaussianly  distributed  displacement  tu(o) .   It 
follows  that 

•3/2 


(6.5)    G(x,t) 


4rt(Kt  +  I  VqV)     expl -1x|  /(hKt  +   2Vq  t^) 


(t  «  l^/v^), 


If  \lf(x,t)  represented  temperature  fluctuation  instead  of  marked  particle 
concentration,  and  k   were  kinematic  heat  conductivity,  then  B  would  be 
called  a  Peclet  number. 
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and  that 

(6.6)      (^/)  =   rx^^G(x,t)d\  =  2Kt  +  y^^t^       (t  «  i^/v^). 


where  x  is  any  component  of  displacement.   We  see  from  (6.5)  and  (6.6) 
that  molecular  diffusion  is  negligible  for  most  times  if  B  »  1,  but 
that  it  dominates  for  very  short  times  t  «  k/v       =  B   ^J^ci' 

For  t  »  ^^/^f^}   ■we  may  anticipate  that  the  convection  acts  as  a 
random  walk  process  with  a  step  length  approximately  I      and  an  effective 
velocity  approximately  v  .    Hence^  we  shall  assume  that  G(x,t)  has  a 


Gaussian  form  in  this  range  also: 


(6.7)     G(x,t)  =  (l+nav^tQt)"5/2exp 


-Ixl'Aav^t^t 


'        (t  »  V^o^' 


(6.8)     (x/)  =  Sav^t^t,  it»l^/y^). 

Here  a  is  a  number  fixed  by  the  shape  of  U. .(x-x'),  and  we  neglect  mole- 
cular  diffusion  completely. 

On  the  basis  of  (6. 5) -(6. 8)  we  shall  distinguish  three  asymptotic 
time  ranges 

Very  short  times:    t  «  B   ^q/^q 

Short  times:   ^''^^q'^'^q   «  ^   «  ^o'^^O 
Long  times:   t  »  ^^A„ 


The  three  ranges  will  be  distinct  if  B  is  large  enough. 


5 ~~ ~~~ 

For  a  static  velocity  field,  the  asymptotic  correspondence  of  random 

convection  to  a  random  walk  does  not  appear  to  hold  in  two-dimensional 

flows . 
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6.2  Iteration  expansion 

Equation  (6.2)  may  be  solved  formally  by  iteration^  in  analogy  to 
the  procedure  discussed  in  Section  5.2.   If  the  result  is  averaged  over 
the  normal  distribution  of  u(x)^  we  find 

t    s 
(6.9)     G(x,t)  =  GQ(x,t)  +   Tds  Tds-  j   d\jd\^G^{x-Y,t-s) 

0    0 


SG_(y-ySs-s')  BG_(y',s') 
XU.  .(y-y-)    °-  ~ 5-^= +  ...    , 

^^  --^  Sy  Sy.  • 


0 


where 


(6.10)    G_(x,t)  =  (ifrt/<t)"^/^exp(-|x|^A/<t) 


is  the  solution  of  the  homogeneous  equation  with  u(x)  =  0.   Equation 
(6.9)  also  may  be  written  as  the  differential  equation 

(6.11)     Q^  -   /c  v2^G(x,t)  =  f  <^  f  <i\\^i^-Y) 


SG  (x-y,t-s)  aG-(y,s) 
X    ^  ~"^ -ii-:= +  ...  ,     G(x,0)  =  5^(x). 

ax^  ay. 

The  higher  terms,  which  are  not  shown  explicitly  in  (6.9)  and  (6.11), 
involve  derivatives  of  U. .  as  well  as  U. .  itself.   If  we  were  to  adopt 

appropriate  characteristic  units,  as  was  done  in  Section  5,  we  would  see 

2 
that  (6.9)  is  an  expansion  in  powers  of  B   (odd  powers  of  B  do  not 

appear) . 
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For  t  «  I     /k=   B   I   /y  y    which  can  include  all  three  of  our  asymp- 
totic time  ranges  if  B  is  large^  G  (x, t)  is  negligible  imless  |x|  «  t    . 

2 
:rx  -L ij^-'    i/- 

explicitly  shown  terms,    and  thereby  obtain 


Then  ve  may  take  U.  .(y-y')  :x   U.  .(O)   =   6    .v        in   (6.9)^    retain  only  the 

X,]  «v   *rt,  X -J  X -J   u 


(6.12)  G(x,t)  =  GQ(x,t)  +  ~   v^^t^  ^\(^'^^  (*  «  ^o^q/^'o^' 

(6.13)  (x/)  =  2Kt  +  v^^t^  (t  «  Bq^q/vq) 

as  the  first  iteration  approximation.    For  t  «  B   ^J^ci'    ^^^   second 
terms  on  the  right  sides  of  (6.12)  and  (6.15)  are  small  corrections,  and 
the  iteration  approximation  is  a  good  one.   For  longer  times,  (6.12) 

provides  a  very  poor  approximation  to  (6.5)  aJ^d-  (6.7)-   It  indicates  that 

1/2 
G(x,t)  is  negligible  for  [x]  »  [Kt)    '    ,  whereas,  according  to  (6.5)  or 

(6.7),  it  should  be  spread  over  the  much  greater  range  |x|  -^  v  t  or 

1/2 
|x|  ~  (av  t  t)  '  .   Moreover,  G(x,t)  is  strongly  negative  in  a  region 

~        1/2  "" 

|x|  <  (Kt)  '   surrounding  the  origin,  in  contradiction  of  the  physical 

requirement  that  the  probability  of  finding  a  particle  be  non-negative. 


It  is  of  Interest,  however,  that  (6.I3)  agrees  with  the  exact  solution 

.  /   6 
throughout  the  short  time  range  t  «  ^qI^q' 

We  wish  now  to  examine  the  convergence  properties  of  the  full  itera- 


tion expansion,  for  B  »  1.   The  expansion  can  be  generated  very  simply 
in  the  short  time  range  if  we  rewrite  (6.5)  in  operational  form  and  then 


To  obtain  (6.12)  we  perform  partial  integrations  and  note  that  G  Is 
its  own  convolution  in  space-time. 

The  fact  that  ^x   )  grows  properly  (t  «  ^J^r.)}    although  G(x,t) 
remains  concentrated  near  the  origin,  is  related  to  the  Indefiniteness 
of  G(x,t).   Although  the  total  probability  is  conserved,  the  total  pro- 
babllity  in  the  region  where  G(x^t)  is  positive  grows  with  t. 
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GQ(x,t) 


expand  in  povers  of  v   .      We  have 

(6.lif)  G(x,t)   =  exp^l  VqV  v2"^GQ(x,t) 


00  . 

n=l       ^ 


{t«l^/.^). 


It  can  be  seen  from  (6.1^)  that  successive  terms  in  the  expansion  display- 
more  and  more  complicated  oscillations  in  sign.   For  B  »  1^  the  succes- 
sive terms  grow  in  magnitude^  except  for  very  short  times.   Thus,  the 
iteration  approximations  to  G(x,t)  do  not  form  a  usefully  convergent 
sequence.  We  should  note,  however,  that  ^x    /is  given  correctly  for 
m  i  N,  in  the  short  time  range,  if  (6.1^+)  is  truncated  after  the  term 
n  =  K. 

The  actual  radius  of  convergence  of  (6.1^)  is  easy  to  determine. 
If  we  examine  the  right  side  of  (6.5)  as  a  function  in  the  complex  v 

plane,  wee  that  it  is  analytic  at  the  origin  but  has  essential  singular- 

1/2 
ities  at  V  =  ±  i(2/c/t)  '  .   Hence,  the  domain  of  convergence  of  (6.1^) 

as  a  power  series  in  v  Is  given  by 

t  <  2k/v/  =  2B^-ll^/v^. 

In  the  case  B  «  1,  which  we  have  not  considered,  this  result  suggests 
that  the  domain  of  convergence  extends  into  the  long  time  range  t  »  ^^^/v  - 
However,  (6.5)  Is  not  valid  there.  Actually  the  iteration  expansion  for 
G(x,t)  appears  to  converge  for  all  t,  if  B  «  1. 
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6.5.  Ciunulant- discard  approximations 

In  analogy  to  Section  5-5^  we  may  constinct  a  sequence  of  approxi- 
mations to  G(x, t)  by  truncating  a  cumulant  expansion  of  the  joint  distri- 
bution of  u(x)  and  \jf(x,t).   This  has  been  carried  out  for  a  more  general 
case  by  Roberts  [^24]  .   The  sequence  of  closures  relevant  to  the  present 
problem  is 

(u  (x')u  (x")^(x,t))  =  U  (x'-x")(,l/(x,t))  , 

(6.15) 

(u.(x')u.(x")ujx'")^(x,t))  =  U.  .(x'-x")  <u^(x"')t(x,t)) 

+  U^Jx'-x"')(u  (£)^(x,t))  +  U   (x"-x"')  (u.(x')t(x,t))  , 


By  using  the  first  of  these  equations  to  close  off  the  pair  of  moment 
equations 


(6.16)    ^Ib  -  '<^^)  (t(x.t))  =  -  /u^(x) 


^x. 
1 


(6.17)    Qt-  ><^T)  {\^i'^^^l'^'>)    =-     \^it')^j 


St(x,t)' 

(x) /   , 

hx. 

J 


we  obtain,  after  manipulation,  the  first  cumulant-dlscard  approximation 


t 
(6.18)     r^  -'^^G(x,t)  =  JdsJd\^\J^.{x-x') 


SG  (x-x%t-s)   dG(x',s) 

X   — "    -—^ ,      G(x,0)  =  6^(x) 

dx.         5x . '  "" 
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It  will  be  noted  that  (6.I8)  bears  a  resemblance  to  the  first  itera- 
tion approximation^  as  given  by  the  explicitly  shown  terms  in  (6.11). 
Equation  (6.l8)^  however,  is  an  homogeneous,  linear,  integro-dlfferential 
equation  instead  of  an  1 nhomo gene ous ,  linear,  differential  equation. 

For  t  «  Bq^qA„j  (6.18)  may  be  reduced  to  the  form 


(6.19; 


(jt   -  "^y  ^^^'^^   =  \^    f    ^^G(x,s)ds,      G(x,0)  =  b^{x), 


of  which  the  solution  is 


(6.20)    G(x,t)  =  exp(/cV  t) 


hnlxi   ai 


&(|^l  -  v^t) 


.  .(4.Kt)-5/2j£t_l 


ax 


exp 


Hi^l-V 


2n 


k-Kt 


+   exp 


(lx|  +  v^t) 


2i 


k-Kt 


For  t  «  B   ^q/^o'    '^^^   range  where  the  first  iteration  approximation  (6.12) 
is  valid,  it  may  be  verified  that  (6.20)  reduces  to  (6.12),  as  would  be 
anticipated  by  argiiments  similar  to  those  in  Section  5.2.   For  longer  times, 
where  the  first  iteration  approximation  is  poor,  (6.20)  also  Is  a  poor 
approximation  to  (6.5)  and.  (6.7).   For  such  times,  (6.20)  gives  to  G(x,t) 

the  form  of  an  outgoing  spherical  wave-front,  consisting  of  a  narrow 

1/2 
positive  crest  [width  ~  (Kt)  '  J  followed  by  a  similar  negative  crest^ 

which  moves  outward  at  a  constant  velocity  v  .   Thus  It  gives  negative 

probabilities,  as  did  the  Iteration  approximation.   The  behavior  of  the 

present  approximation  is  somewhat  better,  however,  in  that  G(x,t)  at  least 
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does  not  remain  concentrated  about  the  origin.   Equation  (6.20)  gives 
for  ^x   )  the  result  (6.15)  obtained  in  the  first  iteration  approxi- 
mation. As  we  have  noted,  this  expression  is  correct  for  t  «  ^„/^„ 
but  not  for  t  »  -f-^^/v  • 

The  next  few  approximations  in  the  cumulant- discard  sequence  are 
not  difficult  to  work  out  in  the  short  time  range  t  «  ^r^/^^'      We  shall 
give  the  final  results  for  the  second  and  third  approximations.   Using 
the  second  of  equations  (6. 15)  to  close  off  the  first  three  equations 
in  the  moment  sequence  which  commences  with  (6.l6)  and  (6.17)^  one  finds 


(6.21)    G(x,t)  =  exp(^V  t) 


6n 


5(|x|) 


a8(|x|  -  V5  V  t) 


12n  X 


alxl 


The  third  cumulant -discard  approximation  is 


(6.22)    G(x,t)  =  -  ^  exp(/c  V^t)  -((5+a/6)|x 


-1 
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■(5-V5)^/V] 


Sixi 


+  (5  -\/S)|x 


_^  d6[|x|-(5+V^)^/\t] 


alxl 


It  is  apparent  that  the  unphysieal  characteristics  of  wave-like  pro- 
pagation and  negative  probabilities  persist  in  these  higher  approximations. 
However,  there  is  a  kind  of  improvement  which  also  characterized  the  higher 
iteration  approximations:   The  quantity  /x    /  ^^  given  correctly  for 
m  S  N  in  the  Wth  cumulant-discard  approximation,  provided  t  «  -^^Aq-   For 
t  »  -t  /v  ,  the  behavior  of  the  second  and  third  cumulant-discard  approxi- 
mations is  complicated,  but,  like  the  first  one,  they  do  not  appear  to  go 
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over  into  the  random  valk  behavior  described  by  (6.7)  and  (6.8).   Instead, 
they  appear  to  yield  ^x   )  oc  t  in  this  time  range. 

The  results  we  have  cited  seem  clearly  to  indicate  that  neither  the 
magnitude  nor  the  importance  of  the  cumulants  of  the  joint  distribution 
of  u  and  \|r  decreases  with  increasing  order  when  B  »  1.  Apart  from  the 
question  of  the  quantitative  accuracy  with  which  the  cumulant -discard 
approximations  represent  G(x,t),  it  is  dlsqiiietlng  that  they  lead  to  nega- 
tive probabilities  for  finding  a  particle. 

Similar  unphysical  results  occur,  and  are  more  disconcerting,  if  we 
use  the  cumulant-discard  approximations  to  trace  the  evolution  of  the 
manifestly  positive-definite  quantity  (  [■i|;(x,t)3  /,  We  shall  outline  this 
analysis  for  the  case  of  pure  convection  (k  =  O).   In  this  case, 

/  [\|/(x, t)]  d  x,  which  plays  a  role  analogous  to  total  kinetic  energy  In 
turbulence  dynamics,  is  conse2rved. 

Let  us  take 

(6.25)         ^(x,t)  =  ^'  +  t'(x,t), 

where  ■ij;  is  a  constant  background  concentration  and  \|f'(x, t)  is  a  fluctua- 
tion which  satisfies  /\j;'(x,t)>  =0   [a  condition  preserved  by  (6.2)3. 
We  wish  to  trace  the  evolution  of 


(6.24)         Y(x,t)  =  ([^•(x,t)]2) 


under  the  condition  that  \|f'(x,0)  is  statistically  independent  of  u. 
From  (6.2),  with  /c  =  0,  we  have 


(6.25)        ^  &'(x,t)]2  =  -u^(x)  J-  [,^'(x,t)]^ 
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a  linear  equation  identical  with  that  obeyed  by  \|f(x,t)  Itself.   The 
sequence  of  eumulant-dlscard  closures  appropriate  to  the  present  prob- 

lems  are  identical  with  (6.I5)  if  we  replace  \(f(x^t)  there  by  [\if '  (x, t)^  . 

7 
Thus  the  first  closure  is 

(6.26)    (u.(x')u  (x")[,|f'(x,t)]^)  =  U,  .(x'-x")Y(x,t). 

It  follows  immediately  from  these  similarities  that  the  cumulant- discard 
approximations  to  f(x^t)  and  G(x^t)  are  related  by 


(6.27)    T(x,t)  =   /  G(x-xSt)^(x',0)d\', 


which  holds  also  if  no  approximation  is  made.   In  particular^  to  obtain 
the  first  cumulant -discard  approximation  for  ¥(x,t),  we  substitute  (6.20)^ 
with  K=  0,  for  G(x,t)  in  (6.27).   This  result  holds  for  all  t^  since 


Bq  =  0O. 


To  illustrate  our  result,  let  us  take 


(6.28)    f(x,0)  =  (2rtd^)'^/^exp 


llxlW 


where  d  is  a  measure  of  initial  spread.   In  the  first  cumulant -discard 
approximation,  we  find 


P 

(6.29)      ^(x,t)  =  -(2rtd2)-5/2  ^ a_ 

"~  Ixl   d  Ixl 


exp  - 


+  exp 


'Moments  of  the  form  ( u^(x' )\|r(x,t))  vanish  if  (i|/'(x,0))  =  0. 
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This  predicts  the  development  of  negative  values  for  l{x,t).      Moreover, 
we  find  from  (6.29)  that,  while  /  f(x,t)d  x  is  conserved  as  in  the  exact 
solution,   /  |Y(x,t)|dx  grows  indefinitely. 

The  closure  (6.26)  is  identical  with  that  employed  in  a  slightly 
more  general  form  by  Reid  [22\    to  investigate  the  evolution  of  the  wave- 
number  spectrum  of  an  isotrQpically  distributed  field  ij/'(x,t)  convected 
by  Isotropic  turbulence.  Reid  did  not  assumfe  a  static  velocity  field, 
as  we  have  done,  but  instead  took  a  field  obeying  the  inviscid  Navier- 
Stokes  equation.   It  is  easily  verified  that  our  result  (6.29)  for  the 
spread  of  an  inhomogeneous  scalar  field  is  unchanged  in  this  case.   The 
introduction  of  viscosity  does  produce  a  change,  but  the  unphysical 
behavior  of  Y(x,t)  persists.   Our  present  results  suggest  that  the  cimiu- 
lant-dlscard  approximations  for  an  isotropically  distributed  scalar 
field  may  lead  to  negative  spectral  densities  for  some  wavenumbers.   The 
explicit  integration  of  Reid's  spectral  equation  would  shed  light  on  this 
question. 

6.^.   Stochastic  model  approximations 

The  first  stochastic  model  of  the  convection  problem  [l6j  has  been 
investigated  by  Roberts  f25j  ■   The  equation  for  G(x, t),  which  we  shall 
not  derive  here,  is 


(6.30)     Q^  -  /<V^^-G(x,t)  =  J  ds  J  d\\J^.{x-y) 


2\     ,        ,  r         r    3  ,    ^   SG(x-y,t-s)  aG(y,s) 

Sx.       by. 
1         J 


G(x,0)  =  6^(x) 


ho 


It  will  "be  noticed  that  (6.50)  resembles  (6.I8)  but  differs  by  being 

nonlinear  in  G(x_,t).   For  t  «  -f-^/v^  U  .(x-y)  in  (6.50)  may  be  replaced 

2 
by  U.  .(0)  =  5  .V   ,   Then  the  equation  may  be  reduced  by  partial  integra- 
tion to  the  form 


(6.51)    0^  -  K^^jGix,t)   =  v^^  J  ds  J  d.\G{x-Y,t-s)^^G{y,s), 

0 

G(x,0)  -  6^(x), 


of  which  the  solution  is 


(6.32)    G(x,t)  =  (2nvQt)"^exp(/<V^t; 


.l(lx|/2vQt)(W^2t2-|^l'')'^^^ 


where 


Ti(s)  =  1     (s  <  1),     Ti(s)  =0    (s  s  1) 


Equation  (6.52)  represents  a  substantial  improvement^  over  any  of 
the  iteration  or  cumulant-discard  approximations  we  have  examined^  in 
that  G(x^t)  is  positive  definite  and  is  spread  out  over  a  region  corre- 
sponding to  the  spread  of  the  exact  solution  (6.5).   In  common  with  the 

/   2\ 

previous  approximations,  the  present  one  gives  {^x       /     correctly  in  the 

range  t  «  I   /v     which  we  are  considering.   Apart  from  these  virtues, 
however,  the  approximation  is  a  poor  one.   Instead  of  falling  off  mono- 
tonically  with  the  increase  of  |x|,  G(x,t)  is  minimum  at  the  origin, 
according  to  (6.32),  and  rises  to  a  cusp  and  cut-off,  at  |x|  =  2v  t, 
which  become  Infinitely  sharp  as  K  -*  0. 

The  behavior  of  the  first  stochastic  model  for  long  times  is  more 
satisfactory,  and  it  represents  a  profound  improvement  over  the  previous 


hi 


approximations.      For  t  »  ^^/"v_.    It   can  be   shown    [25J    that  G(x,t)   does 
take  the   form  (6.7)  •      The  n\:iraber  a  is   given  by 

(6.55)  a  =   (Jv^t^)"^    Jdsjd.\  U.^(y)G(y,s). 

0 

In  order  to  evaluate  this  expression  explicitly,  we  must  integrate  (6.50) 
to  find  G(y, s)  for  all  time  ranges.  The  dominant  contributions  to  (6.53) 
come  from  the  region  ( |y|  ~  I   ,    s  ~  t  /v  ) .   How  accurate  a  quantitative 

u^      \J         u   (J 

approximation  (6.33)  represents  is  not  known. 

The  behavior  of  G(x, t)  has  been  partially  investigated  for  a  tenta- 
tive second  stochastic  model  [16] .   The  approximation  to  the  exact  solu- 
tion is  very  substantially  improved  for  t  «  ^J^r-^-      Outside  of  the 
asymptotic  time  ranges,  the  equations  are  very  complicated. 

6.3.   Formal  relations  among  the  approximations 

We  stated  in  Section  5  that  the  .cumulant -discard  approximations  and 
the  stochastic  model  approach  both  represented  well-defined  summations  of 
selected  terms  from  all  orders  of  the  iteration  expansion.   This  is  also 
true  in  the  present  problem,  and  the  relations^  in  fact,  take  a  simpler 
form.   The  iteration  expansion  of  G(x,t)  consists  of  a  sxm  of  functionals 
of  U.  .(x)  of  ascending  degree.   The  contribution  of  nth  degree  arises  as 
the  average  of  a  product  of  2n  factors  u.   Therefore,  according  to  the 
rules  for  evaluating  moments  of  a  normal  distribution,  this  contribution 
consists  of  (2n)f/2  ni  tenns,  one  for  each  possible  grouping  of  the  u 
factors  into  pairs.   In  the  cumulant -discard  and  stochastic  model  approxima- 
tions, some  of  these  terms  are  discarded  and  others  retained  unaltered. 
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We  shall  not  state  the  selection  rules  here.    Instead,  we  have  indi- 
cated in  Table  1  the  total  number  of  terms  of  each  order  (degree  in  U . . ) 
retained  in  each  approximation.   This  is  a  particularly  interesting  measure 
for  t  «  t  /v  and  k  =   0-      Then  all  the  (2n)l/2  nl  terms  in  the  nth  order 
of  the  iteration  expansion  give  exactly  equal  contributions  to  G(x^t). 


TABLE  1 


Approximation 

1st  Cumulant- discard 
2nd  Cumulant -discard 
5rd  Cumulant -discard 
1st  Stochastic  Model 


Terms  retained  of  order  n 

1 
,n-l 


(2n):/nl(n  +  1)1 


It  may  be  inferred  from  Table  1  that  mere  comprehensiveness  of  summa- 
tion is  not  necessarily  a  good  measure  for  the  validity  of  an  approximation. 
The  third  curaulant-discard  approximation  contains  more  terms  for  large  n  than 
the  first  stochastic  model.   (The  asymptotic  numbers  are  of  order  (5+V^) 
and  n~   '    k^ ,    respectively.)   However,  it  gives  a  poorer  approximation  at  the 
long  times  where  the  high  orders  can  contribute  most  strongly. 

We  should  like  to  mention  here  that  the  structure  of  the  iteration  series 
for  G(x,t),  and  the  infinite  summations  of  terms  from  it  which  our  approxima- 
tions represent,  have  analogs  in  quantum  field  theory.   In  particular,  the 
cumulant-discard  approximations  are  close  analogs  to  the  Tamm-Dancoff  approxi- 
mations for  the  "one-nucleon  propagator"  in  meson  theory  f26j . 


The  iniles  for  the  stochastic  models  are  given  in  [l6]  and  [25J 


h^ 


7 .   CONCLUSION 

We  have  been  concerned  principally  with  a  particular  form  of  the 
closure  problem  of  turbulence  theory:  Given  a  simple  initial  distribution 
for  the  velocity  field  (a  normal  distribution),  how  can  one  obtain  a  sim- 
ple analytical  description  of  the  distribution  at  later  times?  We  have 
indicated  that  a  simple  exact  description  is  not  known  and  is  likely  to 
be  very  difficult  to  find.   This  has  led  us  to  consider  approximation 
approaches  to  the  closure  problem,  and  we  have  examined  several  in  the  con- 
text of  isotropic  turbulence. 

The  most  straightforward  approximation  scheme  we  considered  was  based 
on  expansion  of  the  moments  of  the  distribution  in  powers  of  the  turbulent 
Reynolds  number  R  of  the  initial  state.   This  scheme  can  be  expected  to 
yield  good  approximations  only  for  the  very  restricted  case  R  «  1.  We 
have  pointed  out,  however,  that  the  Reynolds  number_,  or  iteration,  expan- 
sion is  of  central  theoretical  interest  at  all  R  as  a  formal  standard  of 
comparison  for  other  approximations. 

The  second  approach  we  examined  was  the  expansion  of  the  distribution 
in  cumulants.   It  is  more  intimately  and  logically  based  on  the  statistical 
nature  of  the  turbulence  problem.   The  approximations  arising  from  the 
cumulant  expansion  have  not  yet  been  exploited  sufficiently  to  determine 
their  validity.   We  have  seen,  however,  that  the  consequences  of  these 
approximations  in  the  problem  of  convection  by  a  prescribed  random  velocity 
field  do  not  provide  a  favorable  prognosis  for  their  validity  in  turbulence 
dynamics.   The  unfavorable  results  include  negative  probabilities  and  nega- 
tive mean-square  field- amplitude  fluctuations.   They  suggest  that  the 
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cumulants  of  the  velocity  distrilDution  in  turbulence  dynamics  may  not 
decrease  in  importance  with  increasing  order,  when  R  is  large.   It  is 
very  possible,  however,  that  the  convection  example  is  misleading. 
Although  it  presents  a  closure  problem  similar  to  that  of  turbulence 
dynamics,  there  is  the  essential  difference  that  the  equations  are  linear 
in  the  dynamic  variables.   Whether  nonlinearity  in  the  dynamics  serves  to 
reduce  or  to  accentuate  the  unphysical  consequences  of  truncating  the 
cumulant  expansion  remains  to  be  determined. 

If  it  turns  out  that  the  cumulant  expansion  does  not  provide  satis- 
factory closure  approximations  for  isotropic  turbulence  at  high  R  ,  the 
outlook  for  finding  acceptable  alternative  approximations  of  equal  sim- 
plicity does  not  seem  very  favorable.   This  is  illustrated  by  the  third 
scheme  we  have  considered,  that  of  solving  model  problems  which  lead  to 
exact  closed  moment  equations.   Negative  probabilities  and  related 
troubles  cannot  arise  in  this  scheme,  but  the  price  paid  for  the  assurance 
is  a  quite  substantial  increase  in  the  complexity  of  the  final  equations 
which  must  be  solved. 

In  both  the  cumulant  expansion  scheme  and  the  model  scheme,  the 
complexity  of  the  approximations  increases  so  rapidly  with  order  that  only 
the  lowest  approximations  give  much  hope  of  practical  utility.   If  they 
prove  inadequate,  the  need  for  a  genuinely  fresh  theoretical  approach 
would  seem  strong.   Perhaps  there  is  an  alternative  representation  of  the 
statistical  problem  which  leads  to  simpler  and  more  natural  approximations 
than  those  based  on  moments .   Unfortunately,  there  seem  to  be  only  meager 
clues  in  the  branches  of  statistical  physics  which  have  yielded  more 
graciously  to  analytical  treatment.   Turbulence  at  high  Reynolds  numbers 
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is  characterized  by  both  strong  nonlinearity  and  strong  dissipation. 
This  appears  to  be  a  particularly  intractable  combination  of  properties. 

We  should  like  to  close  with  the  remark,  made  earlier,  that  iso- 
tropic turbulence,  despite  Its  geometrical  simplicity,  may  constitute  a 
particularly  severe  testing  ground  for  closure  approximations.  As  we 
have  seen,  the  closure  problem  arises  solely  from  the  self -interaction 
of  the  fluctuating  part  of  the  total  velocity  field.   In  shear  flow  prob- 
lems, the  interaction  of  mean  velocity  field  and  fluctuating  field  makes 
a  major  contribution  to  the  total  dynamics.   Perhaps  this  situation  is 
one  in  which  cruder  closure  approximations  are  permissible.   None  of  the 
approximation  schemes  we  have  discussed  is  limited  to  homogeneous  turbu- 
lence, and  it  is  highly  desirable  that  they  be  investigated  for  more 
general  flows . 
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